Dynamics of Majority Rule in Two-State Interacting Spins Systems 
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We introduce a 2-state opinion dynamics model where agents evolve by majority rule. In each 
update, a group of agents is specified whose members then all adopt the local majority state. In the 
mean-field limit, where a group consists of randomly- selected agents, consensus is reached in a time 
that scales IniV, where N is the number of agents. On finite-dimensional lattices, where a group 
is a contiguous cluster, the consensus time fluctuates strongly between realizations and grows as a 
dimension-dependent power of N. The upper critical dimension appears to be larger than 4. The 
final opinion always equals that of the initial majority except in one dimension. 
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In this letter, we introduce an exceedingly simple opin- 
ion dynamics model - majority rule (MR) - that exhibits 
rich dynamical behavior. The model consists of N agents, 
each of which can assume the opinion (equivalently spin) 
states +1 or —1, that evolve as follows: 

• Pick a group of G spins from the system (with G an 
odd number). This group could be any G spins, in 
the mean-field limit, or it is contiguous, for finite- 
dimensional systems. 

• The spins in the group all adopt the state of the 
local majority. 

These two steps are repeated until the system reaches a 
final state of consensus. While the MR model ignores 
psycho-sociological aspects of real opinion formation pj , 
this simple decision-making process leads to rich collec- 
tive behavior. We seek to understand two basic issues: 
(i) What is the time needed to reach consensus as a func- 
tion of N and of the initial densities of plus and minus 
spins? (ii) What is the probability of reaching a given 
final state as a function of the initial spin densities? 

To set the stage for our results, we recall the corre- 
sponding behavior in the classical 2-state voter model 
( VM) [J , where a spin is selected at random and it adopts 
the opinion of a randomly-chosen neighbor. For a system 
of N spins in d dimensions, the time to reach consensus 
scales as N for d > 2, as N In N for d = 2 (the critical 
dimension of the VM), and as N 2 in d = 1 Because 
the average magnetization is conserved, the probability 
that the system eventually ends with all plus spins equals 
the initial density of plus spins in all spatial dimensions. 

The MR model has the same degree of simplicity as 
the VM, but exhibits very different behavior. Part of the 
reason for this difference is that MR does not conserve 
the average magnetization. Another distinguishing trait 
of MR is the many-body nature of the interaction. This 
feature also arises, for example, in the Sznajd model Q, 



where two neighboring agents that agree can influence a 
larger local neighborhood, or in Galam's rumor forma- 
tion model ||, where an entire population is partitioned 
into disjoint groups that each reach their own consensus. 
The updating of an extended group of spins was also con- 
sidered by Newman and Stein [(| in the Ising model with 
zero-temperature Glauber kinetics 0]. 

We now outline basic features of the MR model. In 
the mean-field limit we give an exact solution for the 
approach to consensus, while for finite dimensions we give 
numerical and qualitative results. 

Mean-field limit. Consider the simplest case where ar- 
bitrary groups of size G = 3 are selected and updated at 
each step. To determine the ultimate fate of the system, 
let E n denote the "exit probability" that the system ends 
with all spins plus when starting with n plus spins. Now 
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is the probability that a group of size 3 has j plus and 
3 — j minus spins in an A-spin system that contains n 
plus spins. The group becomes all plus for j = 2, it 
becomes all minus for j ' = 1, while for j = or 3 there is 
no evolution. Thus E n obeys the master equation [|j 
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which simplifies to 

(n - l)(£„+i - E n ) = (N -n - l){E n - £ n _ x ). (2) 

Writing D n = E n+ i — E n , Eq. becomes a first-order 
recursion whose solution is 



D n = 



D 



r(n) T(N — n — 1) 



(3) 
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To compute the constant B we use the fact that 
Za<„<2V_2 D n = Sjv-i — Ex = 1, due to the bound- 
ary conditions E\ = and -E/v_i = 1. Thus we find 
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The probability to end with all spins minus is simply 
i?jv_„. Since consensus is the ultimate fate of the system 
E n + Eff—n = 1. 

While the minority may win in a finite system, the 
probability for this event quickly vanishes as N — > oo . In 
the continuum limit n, N — > oo with p = n/A~ < 1/2 the 
exit probability is exponentially small: E n cx X N , with 
X = l/[2p p {\ - p) 1_p ]. Only near n = AT/ 2 does the 
exit probability rapidly increase. Employing Stirling's 
approximation we may recast Q into 



E n -> £{y) 



1 



(5) 



where y = (2n- N)/VN. 

We now study the mean time T n to reach consensus 
(either all plus or all minus) when the initial state con- 
sists of n plus and N — n minus spins. (The time to 
reach a specified final state can also be analyzed within 
this framework.) Similar to the reasoning for the exit 
probability, the equation for T n is Q 
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subject to the boundary conditions To = TV = 0. The 
natural choice for the time interval between elementary 
steps is ST — 3/N, so that each spin is updated once per 
unit time, on average. 

The master equation for U n = T n +\ — T„ simplifies to 

(n - l)U n = (N-n- l)t/„-i - (N ~ * )(JV 7 2) , (7) 

n(N — n) 

with the boundary conditions Uq — 1 and Un—i = — 1- 
Apart from the inhomogeneous term, Eq. Q is identical 
to (2). Thus, we seek a solution in a form similar to J^J: 



U„ 



(N -1)(N -2) 



r(n) T(N — n — 1) 
This transforms Q into the difference equation 
r(n- l)r(JV-n- 1) 



(8) 



(9) 



The symmetry of Eqs. (JJJ and ansatz iJSJ under the trans- 
form n — » N — 1 — n, and the antisymmetry of the 
boundary conditions Uq = 1 and Un-i = — 1 imply that 
= -t/jv-i-n and V„ = -Vw-i-n- 
For concreteness, we shall take A/" to be odd and define 
k = (N — l)/2. Then the above boundary conditions on 




FIG. 1: Consensus time T n versus p = n/N for N = 81, 401, 
2001, 10001, and 50001 (gradual steepening). The curves are 
symmetric about p — 1/2. 



U n and V n imply that V4 = 0. Starting from this value 
and using Eq. © we recursively obtain, for all j < k 
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r(fe - i) Y{k + i - 1) 
- (fc - i + 1) (fc + i) 



(10) 



For n > 1, the average time T n = Ti + ^7 becomes 



T n = l + 2k(2k- 1)^2 



lT(j)T{2k-j) 



(11) 



with Vj given by Q1U[I. 

For the maximal time T max = (with fc = (A 7 " — 1) /2), 
we obtain 



T max = 1 + 2fc(2fc - 1) ^ 5, 



k.rn 



(12) 



with 



Sk,m — ^ ^ 



r(fc + j — m) Y(k + m — j — 1) 



^ r(j) T(2fc — j) (k + j — m + 1) (fc + m - j) ' 

A detailed asymptotic analysis shows that 

T max ->21nAT N-^oo. (13) 

For biased initial conditions, it is convenient to con- 
sider the limits n, N — > 00, but with p = n/N kept fixed 
and distinct from 1/2. Now the leading behavior of T n is 



T n -> 2jfe(2jfe-l)^ 



1 



- (2*c-j-l)(i + l)(2fc-j) 



In AT as N->oo. 



(14) 



3 



Thus for biased initial conditions, the consensus time also 
scales as IniV, but with amplitude equal to one. A de- 
tailed analysis |j| indicates that the amplitude sharply 
changes from 2 to 1 within the layer iV -1 / 2 -C p— 1/2 <C 1 
as the system is moved away from the symmetric initial 
condition p = 1/2 (Fig. El) 

One dimension. To implement the dynamics in one 
dimension, we define a group to be G consecutive spins. 
If there is no consensus in the selected group, the opinion 
of the minority-opinion agents are changed so that local 
consensus obtains. We parameterize the opinions by the 
spin states S = ±1. For the simplest case of group size 
G = 3, let S, S', S" be the spins in the group that is being 
updated. Focusing on spin S, this spin flips with rate 



W(S^-S) = (l + S'S") 



1-S 



(S' + S") 



(15) 



The factor l+S'S" ensures that spin S can flip only when 
S' = S" , while the quantity within the square brackets 
ensures consensus after spin S flips. Since S 2 = 1, this 
rate can be simplified to W = 1 + S'S" - S(S' + S"). 

In one dimension, each spin Sj belongs to three groups: 
(Sj-2, Sj-i, Sj), (Sj-i, Sj, Sj+i), and (Sj, Sj+i, Sj+2). 
Therefore the total spin-flip rate is given by 

W(Sj — > — Sj) = 3+jSj_2<Sj— 1 + Sj—iSj+i + >Sj_|_iiSj_|_2 
— Sj (Sj—2 +2Sj—i + 1Sj+i +Sjj r 2)j 

which depends on the state of the two nearest neighbors 
and the two next-nearest neighbors of Sj. The equation 
of motion for the mean spin is 



dt 



-2(SjW(Sj 



(16) 



With the flip rate given above, first-order terms (Sj) are 
coupled to third-order terms (Sj-iSjSj+i). For a spa- 
tially homogeneous system this gives 



dm 1 
dt 



6(mi 



m 3 ) 



(17) 



where to 1 = (Sj) and TO3 = (Sj-iSjSj+i). This coupling 
of different-order correlators makes analytical progress 
challenging. In contrast, the different-order correlators 
decouple in the VM and dmx/dt = 0jQ|- In the mean- 
field limit of MR, TO3 = m\ and the resulting equation 
reproduces the consensus time growing as \rtN and the 
fact that the initial majority determines the final state. 

Despite this distinction with the VM in one dimension, 
the dynamics of the MR can be usefully reformulated in 
terms of the domain walls between neighboring opposite 
spins. As long as walls are separated by at least two 
sites, each undergoes a symmetric random walk, exactly 
as in the VM. However, when two domain walls occupy 
neighboring bonds, then the local spin configuration is 

... 1 ... and these two walls are doomed 

to annihilate. Because of this close correspondence with 
the VM, we expect, and verified numerically, that the 
density of domain walls N(t) decays as t~ x / 2 . 




FIG. 2: Probability that a one-dimensional system of length 
10 4 ends with all spins plus as a function of the initial density 
p of plus spins. Each data point is based on 2000 realizations. 
The dashed line is the corresponding VM result. 



More quantitatively, we study the densities of plus and 
minus domains of length n, P n and Q n , respectively |lfj|. 
The number densities of plus and minus domains are 
identical, N(i) =^P n — J2Qn, while the fractions of 
positive and negative spins are given by L + (t) = ^nf„ 
and L—(t) = ^2,nQ n , respectively (with L + + L_ = 1). 
The equations of motion for these moments are 



dN 
~dt 



-3(Pi + Qi), and 



dL + 
~df 



3(Qi-Pi) (18) 



Substituting N - t^ 1 / 2 into JTSJ gives P x - Q x - t^/ 2 . 
Therefore L+(t)-L+ (00) - i" 1 / 2 . Thus even though the 
fractions of plus and minus spins vary with time, they ul- 
timately saturate to finite values. Correspondingly, the 
exit probability has a non-trivial dependence on the ini- 
tial magnetization in the thermodynamic limit (Fig. 
This should be compared with the VM result E(p) = p 
that follows from the conserved VM magnetization. 

Higher Dimensions. There are many natural ways to 
implement MR in dimension d > 1. One possibility is 
to update groups of three spins at the corners of elemen- 
tary plaquettes on the two-dimensional triangular lattice. 
This gives to a spin-flip rate of a similar form to that in 
Eq. I|15fl and again leads to the equation of motion for 
the mean magnetization being coupled to a third-order 
correlator. On hypercubic lattices, a natural definition 
for a group is a spin plus its 2d nearest-neighbors in the 
d coordinate directions (von Neumann neighborhood). 

On these finite-dimensional lattices, MR differs from 
both the VM and the Ising model with zero-temperature 
Glauber kinetics (IG) 0- On the triangular lattice, for 
example, an elementary plaquette of three plus spins in 
a sea of minus spins cannot grow in the IG model, but it 
can grow in the MR model. Additionally, straight inter- 
faces are not stable in the MR model, but they are stable 
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FIG. 3: Most probable consensus time T mp versus N for zero 
initial magnetization in one dimension (o), on the triangu- 
lar lattice (A), and hypercubic d-dimensional lattice with a 
(2d+l)-site neighborhood for d = 2 (□), d = 3 (*), and d = 4 
(+). The lines represent the best power-law fits with respec- 
tive slopes 2.06, 1.23, 1.24, 0.72, and 0.56. 

in the IG model. On the other hand there is a (small) 
surface tension in the MR dynamics that smooths con- 
vex corners, just as in the IG model. Thus both the MR 
and the IG dynamics lead to relatively compact clusters 
while the VM naturally gives ramified clusters. 

We simulated the time T until consensus is reached 
in one through four dimensions when the system starts 
from zero magnetization. In one dimension, there is only 
a single characteristic time scale, the mean time T, that 
grows as N 2 . The distribution pjv (t) of the scaled time 
t = T/T approaches a well-defined limiting distribution 
Poo(r) for N — ► oo. At late stages of the evolution such 
that only two domain walls remain, we estimate analyti- 
cally and verify numerically that this distribution has an 
exponential long-time tail Poo(t) ~ e~ T as r — > oo. 

In two dimensions (d = 2), we find evidence of at least 
two characteristic times, corresponding to two different 



routes to consensus. In most realizations, one opinion 
quickly becomes dominant and eventually wins. In the 
remaining realizations, however, the system reaches a 
configuration in which opinions segregate into two (or 
very rarely four or more) nearly straight stripes. For 
the IG model, the interfaces between such stripes would 
eventually become perfectly flat and this would be the 
final state 11. However, such stripes are ultimately un- 
stable in the MR model so that consensus is eventually 
reached, albeit very slowly. 

In three and higher dimensions, there are apparently 
numerous characteristic times that all scale differently 
with TV. A similar behavior, associated with a vast num- 
ber of metastable states, was previously observed in the 
three-dimensional IG model [llj . The broad distribution 
of consensus times suggests that the most probable time 
T mp is an appropriate characteristic scale. This quan- 
tity indeed has much more convincing power-law behav- 
ior than the mean time T for d = 2 and 3. We obtain 
Tmp ~ N z with z = 1.23, 0.72, and 0.56 in d = 2, 3, and 
4, respectively, with large fluctuations occurring in the 
most probable consensus time in 4 dimensions (Fig. yj. 

In summary, majority rule leads to rich dynamics. In 
the mean- field limit, there is ultimate consensus in the 
state of the initial majority. The mean consensus time T 
scales as In N, where N is the number of agents in the 
system, and the amplitude of T changes rapidly between 
1 and 2 as a function of the initial magnetization. One 
dimension is the only case where the minority can ulti- 
mately win. Here, single-opinion domains grow as t 1 / 2 
so that T oc A 2 . Because the magnetization is not con- 
served, the probability to reach a given final state has a 
non-trivial initial state dependence. For d > 1, the initial 
majority determines the final state. The consensus time 
grows as a power law in N with a dimension-dependent 
exponent. Large sample-to-sample fluctuations in T arise 
whose magnitude increases with dimension. Mean-field 
behavior is not reproduced in four dimensions, so that 
the upper critical dimension is at least greater than four. 

We thank R. Durrett, T. Liggett, C. Newman, and 
D. Stein for advice and NSF grants DMR9978902 and 
DMR0227670 for financial support of this research. 



[1] See e.g., R. D. Friedman and M. Friedman, The Tyranny 
of the Status Quo (Harcourt Brace Company, San Diego, 
1984); W. Weidlich, Sociodynamics: A Systematic Ap- 
proach to Mathematical Modelling in the Social Sciences 
(Taylor & Francis, London, 2000). 

[2] T. M. Liggett, Interacting Particle Systems (Springer- 
Verlag, New York, 1985). 

[3] P. L. Krapivsky, Phys. Rev. A 45, 1067 (1992). 

[4] K. Sznajd-Weron and J. Sznajd, Int. J. Mod. Phys. C 
11, 1157 (2000); D. Stauffer, J. Artif. Soc. Soc. Simul. 5, 
no. 1 (2002). 

[5] S. Galam, Physica (Amsterdam) 274, 132 (1999); Eur. 



Phys. J. B 24, 403 (2002); cond-mat/ 0211571. 
[6] C. M. Newman and D. L. Stein, Phys. Rev. E 60, 5244 
(1999). 

[7] R. J. Glauber, J. Math. Phys. 4, 294 (1963). 

[8] S. Redner, A Guide to First-Passage Processes (Cam- 
bridge University Press, 2001). 

[9] P. L. Krapivsky and S. Redner, in preparation. 
[10] P. L. Krapivsky and E. Ben-Nairn, Phys. Rev. E 56, 3788 
(1997). 

[11] V. Spirin, P. L. Krapivsky, and S. Redner, Phys. Rev. E 
65, 016119 (2002). 



